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ABSTRACT 

Using the Spitzer IRAC camera we have obtained mid-IR photometry of the 
red giant branch stars in the Galactic globular cluster 47 Tuc. About 100 stars 
show an excess of mid-infrared light above that expected from their photospheric 
emission. This is plausibly due to dust formation in mass flowing from these 
stars. This mass loss extends down to the level of the horizontal branch and 
increases with luminosity. The mass loss is episodic, occurring in only a fraction 
of stars at a given luminosity. Using a simple model and our observations we 
derive mass loss rates for these stars. Finally, we obtain the first empirical mass 
loss formula calibrated with observations of Population II stars. The dependence 
on luminosity of our mass loss rate is considerably shallower than the widely used 
Reimers Law. The results presented here are the first from our Spitzer survey of 
a carefully chosen sample of 17 Galactic Globular Clusters, spanning the entire 
metallicity range from about one hundredth up to almost solar. 

Subject headings: stars: Population II, mass loss - circumstellar matter - in- 
frared: stars 
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Introduction 



Though the current generation of theoretical models can reasonably reproduce the gen- 
eral framework of stellar evolution, there are still a number of physical phenomena which are 
poorly understood. Among these the stellar mass loss (hereafter ML) is one of the most vex- 
ing. This is particularly the case for cool stars. Indirect evidence shows that ML strongly 
affects all of their late stages of evolution, yet we have little theoretical or observational 
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borate ML into our models. We have relied on empirical laws like 
bj) based on observations of Population I giants. While subsequent 



work (lMullanlll978l : lGoldberdll979l : I Judge! 1 1 99 ll : ICatelanl l2000h led to slight refinements, a 



"Reimers Law" or some variant has been the only basis for stellar evolutionary models of 
cool stars at all ages and metallicities. Indeed, a ML law directly calibrated on Popula- 
tion II low-mass giants has never been determined. Galactic globular clusters (GGCs) are 
particularly attractive observational targets because there is so much indirect, but quanti- 
tative, evidence for ML. This includes the observed morphology of the horizontal branch 
(HB) in the cluster color- magnitude diagrams (CMDs), the pulsational properties of the RR 



giant branch (RGB) tip, and th e masses inferred for white dwarfs in GCCs (IRood 
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Direct evidence for ML in GGC RGB stars has been inferred from b l ue sh ifted features 
of photospheric lines like Ha (jDupred Il986l ; iMauas. Cacciari. Pasquinil 120061 ) . These are 



difficult observations and hard to convert to ML rates. One can also observe the outflowing 
gas much further from the star after it has formed dust. We detected this circumstellar (CS) 
matter using ISO C AM and were able to derive ML rates for a modest sample of RGB stars 
(jOriglia et al.ll2002l ). Mid-IR observations have the advantage of sampling an outflowing gas 
fairly far from the star (typically, tens/hundreds stellar radii). Such gas typically left the star 
a few decades previously, hence the inferred ML rate is also smoothed over such a timescale. 
From a combined physical and statistical analysis, our ISOCAM study provided ML rates and 
frequency. We found that i) the largest ML occurs near the very end of the RGB evolutionary 
stage and is episodic, ii) typical rates are in the range 10~ 7 < dM/dt < 10~ 6 Af Q yr -1 , Hi) 
the modulation timescales must be greater than a few decades and less than a million years, 
and iv) there is evidence for dusty shells at even the lowest metallicities. However, our 
ISOCAM survey suffered from two significant limitations. The small sample of observed 
clusters and the consequent low number of detected giants allowed us to reach only weak 
conclusions on the ML dependence on metallicity and HB morphology. Further, the modest 
spatial resolution and sensitivity compromised our ability to measure lower ML rates near 
the RGB tip and made it impossible to explore ML much below the RGB tip. 



Fig. 1. — Three color (3.6/im (blue), 6/im (green), 8/mi (red) mosaiced image of 47 Tuc from 
IRAC. 
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Fig. 2. — Absolute CMDs in the near IR and Spitzer planes of 47 Tuc. Left panel: stars 
detected in the shallow survey down to K « 11 or Mbd ~ over a useful 8' x 5' field of 
view. Right panel: stars detected in the deeper survey down to K ss 14 which excludes the 
inner « 2' (by radius) region. The dashed horizontal line marks the position of the HB level. 
Stars with (K — 8)o color excess at a 3cr level are marked as filled triangles. 
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2. The Spitzer IRAC survey 



A sample of 17 massive GGCs, 4-5 per each 0.5 dex bin in metallicity between [Fe/H] = 
—2.3 and —0.5 and different HB morphologies within each bin has been observed with IRAC 
onboard Spitzer with 26 hr of observing time allocated to our program (ID #20298) in 
Cycle 2. For all these clusters complementary near-IR and UV photometry are available to 
properly characterize both the red and the blue sequences. 

In this Letter we discuss the first cluster analyzed, 47 Tuc. This cluster was observed 
on September 21, 2005. We used a 1 x 3 grid with small cycling dither pattern to cover 
the central 9' x 5' of the cluster in all the four IRAC filters and the High Dynamic Range 
Readout Mode, to avoid saturation of the brightest giants. A total observing time of 1.2 hr 
allowed us to reach K < 14 with S/N ~ 20. Fig. [1] shows the mosaiced three-color image of 
47 Tuc. The Post BCD mosaic frames from the Spitzer Pipeline (Software Vers ion: S13.2.0) 
have been photometrically reduced with ROMAFOT (IBuonanno et al.l Il983l ). a software 
package optimized for Point Spread Function fitting in crowded and undersampled stellar 
fields. Because of the intrinsic luminosity of 47 Tuc, the crowding in its central region is 
critical even at 8 /im. We thus obtained two photometric catalogs: a shallow one down to 
K ~ 11 over a useful 8' x 5' field of view and a deeper one down to K » 14 which excludes 
the inner m 2' (by radius) region. The overall photometric uncertainty in all the four IRAC 
filters is < 0.1 mag. Complementary ground-based near-IR photometry of the central region 
at high spatial resolution has been obt ained using IRAC2 ES02.2m, and SOFI ESO-NTT 
(IFerraro et al.ll2000i ; IValenti et al.ll2004l ) and supplemented with 2MASS data in the external 
region. The degree of completeness of the near-IR survey is 100% over the full magnitude 
range covered by the Spitzer survey. The shallow Spitzer catalog counts almost 400 stars 
and it is > 86% complete in the upper two RGB magnitudes, and ~ 66% complete at 
> Mboi > —2, while the deeper Spitzer catalog in the outer region counts almost 500 stars 
and it is >83% complete down to K < 14. 

The two Spitzer catalogs have been combined and placed onto the 2MASS astrometric 
system by cross-correlating the Spitzer and the ground-based near-IR surveys. The final 
catalog contains almost 800 stars with J, K, and Spitzer photometry in each of the four 
filters. 

The dereddened K magnitudes and (J — K) colors have been also used to compute 
th e stellar bolometri c magn itudes and effective temperatures, by using the transformation 
of iMontegriffo et al.l (119981). reddening [E(B — V) = 0.04) and distance modulus ((m — 
M)o = 13.32) from lFerraro et al.l (119991 ). Fig.[2]shows the absolute color-magnitude diagrams 
(CMDs) of 47 Tuc in the near and mid IR planes. By knowing the stellar temperature and 
bolometric magnitude and using a typical RGB stellar mass of 0.8-0.9 M , we finally estimate 



- 5 - 



the stellar radius and gravity from the following equations: 

U = ^Rl x aT e 4 ff g = GMjRl 

3. Dust excess and mass loss rates 

As CS dust condenses in an outflowing wind, it can be detected as a mid-IR excess. 
In our pilot project using ISOCAM photometry in the 10 /im window and assuming a vB v 
emissivity, we showed that the bulk of CS dust around the RGB tip giants typically has 
temperatures in the range 300-500 K. The IRAC bands between 3.6 and 8 /im can also be 
used to detect this warm dust when coupled with near IR photometry used to properly 
characterize the stellar counterpart. 

In order to select candidate stars with dust excess, we define first the mean ridge lines 
in each of the K , (K — IRAC)® CMDs to define the average color of the stars with purely 
photospheric emission and to determine the overall photometric errors (a) at different mag- 
nitude bins. These (K — IRAC)q colors are ~ 0.0 ±0.1 along the entire RGB range sampled 
by our survey. This is in good agreement with the prediction of theoretical model atmo- 
spheres with T e ff = 3500-5000 K from the Kurucz database. Since the 8 /xm IRAC band is 
the most sensitive to warm dust emission, stars are flagged as dusty if they show a (K — 8)o 
color excess > no . These stars are also the reddest in the other IRAC bands. For these 
stars, we also compute the pure dust emission at 8 /im, by subtracting from the the total 
observed flux the photospheric contribution, given by _pP hot = i^ ega x 10 J -0 - 4 *^ - ^ -8 ^ 110 *)], 
where (K — 8) p hot is the mean ridge line color without dust excess. 

The number of giants with dust excess in the shallow survey of 47 Tuc is 93 at a 3a level. 
In the deeper survey of the cluster outer region, no stars with dust excess have been detected 
below the HB level. The seven 47 Tuc stars which showed dust excess in our ISOCAM survey 
have been also detected by Spitzer and confirmed as dusty stars. A mong the Spitzer dust y 



giants, four are known long period variables (VI, V4, V6, fc V8 in I Clement et al.l (120011 )) 



and an other 17 stars are classified as AGB stars by iBeccari et al.l ( 120061 ). Hence, in the 



following we classify the remaining 74 giants as true RGB dusty stars. 

In order to obtain the ML rates we use our modified version of the DUSTY code 



( llvezic. Nenkova fc Elitzurlll999l ; lElitzur fc Ivezidl200ll ). to compute the emerging spectrum 



and dust emission at the IRAC wavelengths. We adopt Kurucz model atmospheres for the 
heating source and for the dust a mixture of warm silicates with an average grain radius 
a = 0.1 /im. Slightly different choices for the latter two parameters have negligible impact in 
the resulting IRAC colors and mass loss rates. While radiation pressure acting on dust might 
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Fig. 3. — ML rates for the Spitzer sources with dust function of the observed 

(K — 8)0 color (top panel), bolometric magnitude (middle panel) and normalized stellar 
luminosity (bottom panel). In the top panel stars with Mbd < —2 and Mbd > —2 are 
plotted as filled and open squares, respectively. In the middle and lower panels filled big 
triangles mark those giants which are known long period variables, filled small triangles are 
other AGB stars. In the middle panel the evolutionary timescale At in 2 luminosity intervals 
are also reported. In the bottom panel all of the AGB stars are excluded from the fit (solid 
lin es). Typical ran dom error bars are shown in the botto m right corner. The empirical laws 
by Reimers ( 1975al Jbh with 77 = 0.3 (short-dashed line) and Catelan ( 2000 ) (long-dashed line) 
calibrated on Population I giants are also shown for comparison. 
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plausibly drive winds in luminous, metal rich red giants (jWillsoru 120001 ) . the GGC stars are 
generally neither luminous nor metal rich enough for this mechanism to be efficient. Hence 
we run the DUSTY code under the general assumption of an expanding envelope at constant 
velocity v exp with a density profile rj oc r -2 , a dust temperature for the inner boundary r in 
of 1000 K and a shell outer boundary r out = 1000 x r in . We then computed a large grid of 
DUSTY models with stellar temperatures in the 3500-5000 K range and optical depths at 8 
jum (Yg) between 10~ 5 and 10 _1 . For each star with dust excess, we enter the grid with its 
empirical stellar temperature and (K — IRAC)q colors, and we exit with prediction for the 
optical depth, emerging flux, dust fractional contribution and envelope radius. The mass 
loss rates are computed by using the formula 



where 



dM/dt = 47rr 2 ut x p dust x v cxp x 5 



Pdust oc p g r 8 F 8 (obs)/F 8 (mod)L' 2 



is the dust density, p g = 3gcm~ 3 is the grain density, F 8 (obs) and Fg(mod) are the observed 
and model dust emission at 8 micron, respectively, D the distance and 5 the gas to dust 
ratio. A lower limit to 5 is given by 1/Z where Z is the global metallicity. v cxp is a free 
parameter, which should scale like 5 -0 ' 5 if dust and gas are coupled. If the number of grains 
is increased (by decreasing 5), the mo mentum per grain is shared with fewer gas mole cules, 
thus implying an enhancement of v exp (IHabing. Tignon fc Tielenslll994l ; lvan Loonll2000l ) . For 



47 Tuc we adopt 5 ~ 200 and t> exp = 10 kms 1 . The latter is the aver age expansion 



velocity measured in luminous, low mass giants (see e.g. lNetzer fc Elitzurlll993l ) which ranges 
between a few and ~ 20 kms -1 . 



4. Results and Discussion 

For the dusty stars in 47 Tuc sampled by our Spitzer survey, Fig. [3] shows the inferred 
mass loss rates as a function of i) the observed (K — 8)o color, ii) the bolometric magnitude 
and in) the normalized luminosity. The mass loss rate increases with increasing color excess 
and stellar luminosity. Also, the bulk of the mass loss along the RGB should occur above 
the HB level. The provisional empirical law based on such a first set of observations gives: 

dM/dt = C x 4 x 10- 10 x {L/gR)% 4 

where 

C = (5/200) - 5 x (v cxp /10) x (p g /3) 
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and L Q , g Q , and R Q are the stellar luminosity, gravity and radius in solar units. In this 
study of 47 Tuc C — 1. Only true RGB stars are used to derive the fitting formula. Errors 
are as follows: ~ 10% for L/gR and the fit exponent and ~ 25% for the fit zero point. The 
latter is also a good estimate of the average random error on the final mass loss rates. As 
shown in Fig. [3j this new law calibrated on Population II RGB stars is signifi cantly flatter 



than the original Reimers formulation and the one revised by ICatelanl (120001 ). which have 
slopes of 1 and 1.4, respectively. 

As we found in our ISOCAM survey, only a fraction of stars along the RGB are currently 
losing mass: ~ 32% in the upper 2 mag, ~ 16% down to the HB. For this reason we 
conclude that ML is episodic. Basically mass loss is "on" for only some fraction of the time, 
f nn . By using a suitable evolutionary track for a RGB star of M = 0.9 M and Z = 0.004 



( iPietrinferni et al.ll2006l ). we can derive the evolutionary timescale in each luminosity interval 
(see Fig. [3]). Multiplying this by f on , we find that the timescale ML is "on" is less than a 
few Myr in each interval. By using the simple equation 

AM RG b = (dM/dU x At* x / on .) 

to integrate the mass loss formula multiplied by f on over the RGB evolution time At, we 
find that the total mass lost on the RGB is AM RGB ~ 0.23 ± O.O7M . 

In forthcoming papers we will present the results for the other clusters, with the ultimate 
goal of providing the first empirical ML law for Population II stars calibrated over a large 
range of metallicity, and investigating whether observed ML in individual stars within a 
cluster correlates with that cluster's HB morphology and if ML itself is involved in the 
second-parameter problem. 

This research was supported by the Ministero dell'Istruzione, Universita e Ricerca 
(MIUR), the PRIN-INAF 2006 and the Italian Space Agency (ASI). RTR and RMR ac- 
knowledge support from Spitzer Science Center Grant GO-20298. We warmly thanks Elena 
Valenti for having provided the SOFI near IR photometric catalog and the anonymous Ref- 
eree for his/her useful comments. 
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